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Abstract 

A few multiplets that can be added to the SM contain a lightest neu¬ 
tral component which is automatically stable and provides allowed 
DM candidates with a non-standard phenomenology. Thanks to coan¬ 
nihilations, a successful thermal abundance is obtained for well dehned 
DM masses. The best candidate seems to be a SU(2 )l fermion quintu¬ 
plet with mass 4.4 TeV, accompanied by a charged partner 166 MeV 
heavier with life-time 1.8 cm, that manifests at colliders as charged 
tracks disappearing in tt^ with 97.7% branching ratio. The cross sec¬ 
tion for usual NC direct DM detection is dsi = P 1.0- 10“'^^cm^ where 
/ ~ 1 is a nucleon matrix element. We study prospects for CC direct 
detection and for indirect detection. 


1 Introduction 

The Dark Matter (DM) problem calls for new physics beyond the Standard Model (SM). Its 
simplest interpretation consists in assnming that DM is the thermal relic of a new stable 
neutral particle with mass M ~ ~ TeV where Tq ~ 3 K is the present temperature 

of the universe and Gn is the Newton constant. Attempts to address the Higgs mass hierarchy 
problem typically introduce a rich amount of new physics at the weak scale, including DM 
candidates; supersymmetry is widely considered as the most promising proposal [1] . However (i) 
no new physics appeared so far at collider experiments: the simplest solutions to the hierarchy 
problem start needing uncomfortably high fine-tunings of the their unknown parameters [2]; 
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(ii) the presence of a number of unknown parameters (e.g. all sparticle masses) obscures the 
phenomenology of the DM candidates; (ill) the stability of the DM candidates is the result of 
extra features introduced by hand (e.g. matter parity). 

We here explore an opposite, minimalistic approach: focussing on the Dark Matter problem, 
we add to the Standard Model (SM) extra multiplets X + h.c. with minimal spin, isospin and 
hypercharge quantum numbers, and search for the assignments that provide most or all of the 
following properties: 

1. The lightest component is automatically stable on cosmological time-scales. 

2. The only renormalizable interactions of X to other SM particles are of gauge type, such 
that new physics is determined by one new parameter: the tree-level mass M of the 
Minimal Dark Matter (MDM) multiplet. 

3. Quantum corrections generate a mass splitting AM such that the lightest component of 
X is neutral. We compute the value of M for which the thermal relic abundance equals 
the measured DM abundance. 

4. The DM candidate is still allowed by DM searches. 

In section 2 we list the possible candidates. In section 3 we compute the mass splitting. In 
section 4 we compute the thermal relic abundance of X and equate it to the observed DM 
abundance, inferring the DM mass M. In section 5 we discuss signals and constraints from DM 
experiments. In section 6 we discuss collider signals. Section 7 contains our conclusions and a 
summary of the results. 


2 The Minimal DM candidates 


We consider the following extension of the SM: 




■^SM + C 


X{ilp + M)X 
\D^X\‘^ - M‘^\X\^ 


when T” is a spin 1/2 fermionic multiplet 
when T” is a spin 0 bosonic multiplet 


( 1 ) 


where D is the gauge-covariant derivative, c = 1/2 for a real scalar or a Majorana fermion 
and c = 1 for a complex scalar or a Dirac fermion: in all cases we assign X in the minimal 
non-chiral representation of the gauge group, and M is the tree-level mass of the particle. 

We want to identify the cases in which X provides a good DM candidate. Therefore we 
assume the following gauge quantum numbers: 


3. X has no strong interactions [3]. 

2. X is an n-tuplet of the SU(2 )l gauge group, with n = {1, 2, 3,4, 5,...}. 

1. For each value of n there are few hypercharge assignments that make one of the compo¬ 
nents of X neutral, 0 = Q = + Y where T 3 is the usual ‘diagonal’ generator of SU(2 )l. 

For a doublet, n = 2, one needs Y = 1/2. For a triplet, n = 3, one needs Y = 0 (such 
that the component with T 3 = 0 is neutral), or D = 1 (such that the components with 
|r 3 | = 1 are neutral). For a quadruplet, n = 4, Y = {1/2, 3/2}. For a quintuplet, n = 5, 

r = {0,1,2}. 
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Quant 

um numbers 

DM can 

DM mass 

'^DM± ~ '^DM 

Events at LHC 

cTsi in 

SU(2)l 

U(1)k 

Spin 

decay into 

in TeV 

in MeV 

0 

0 

10“^® cm^ 

2 

1/2 

0 

EL 

0.54 ±0.01 

350 

320 d- 510 

0.2 

2 

1/2 

1/2 

EH 

1.1 ±0.03 

341 

160 ^ 330 

0.2 

3 

0 

0 

HH* 

2.0 ±0.05 

166 

0.2 ^1.0 

1.3 

3 

0 

1/2 

LH 

2.4 ±0.06 

166 

0.8 ^4.0 

1.3 

3 

1 

0 

HH, LL 

1.6 ±0.04 

540 

3.0 d-10 

1.7 

3 

1 

1/2 

LH 

1.8 ±0.05 

525 

27^90 

1.7 

4 

1/2 

0 

HHH* 

2.4 ±0.06 

353 

0.10^0.6 

1.6 

4 

1/2 

1/2 

(LHH*) 

2.4 ±0.06 

347 

5.3 d-25 

1.6 

4 

3/2 

0 

HHH 

2.9 ±0.07 

729 

0.01 d-0.10 

7.5 

4 

3/2 

1/2 

(LHH) 

2.6 ±0.07 

712 

1.7 d-9.5 

7.5 

5 

0 

0 

{HHH*H*) 

5.0 ±0.1 

166 

< 1 

12 

5 

0 

1/2 

— 

4.4 ±0.1 

166 

< 1 

12 

7 

0 

0 

— 

8.5 ±0.2 

166 

< 1 

46 


Table 1: Summary of the main properties of Minimal DM candidates. Quantum num¬ 
bers are listed in the first 3 columns; candidates with T 7 ^ 0 are allowed by direct DM searches 
only if appropriate non-minimalities are introduced. The fth column indicates dangerous decay 
modes, that need to be suppressed (see sec. 2 for discussion). The 5th column gives the DM 
mass such that the thermal relic abundance eguals the observed DM abundance (section 4)- The 
6th column gives the loop-induced mass splitting between neutral and charged DM components 
(section 3); for scalar candidates a coupling with the Higgs can give a small extra contribution, 
that we neglect. The 7th column gives the 3a range for the number of events expected at LHC 
(section 6). The last column gives the spin-independent cross section, assuming a sample vale 
/ = 1/3 for the uncertain nuclear matrix elements (section 5). 

For each potentially successful assignment of quantum numbers we list in table 1 the main 
properties of the DM candidates. 

The ‘decay’ column lists the decay modes into SM particles that are allowed by renormaliz- 
ability, using a compact notation. For instance, the scalar doublet in the hrst row can couple as 
XiLji3E°‘eijeai3 where L is a SM lepton doublet, E is the corresponding lepton singlet, i,j are 
SU(2)L-indices, a, (3 are spinor indices, and e is the permutation tensor; therefore the neutral 
component of X can decay as Xq —»• ee. For another instance, the fermion doublet in the second 
row can couple as XaiEgHjEtjEap, where H is the Higgs doublet: its neutral component can 
decay as Xq —>• eh. 

In general, one expects also non-renormalizable couplings suppressed by 1/A^ (where A 
is an unspecihed heavy cut-off scale, possibly related to GUT-scale or Planck-scale physics). 
These give a typical lifetime r ~ A^^ TeV~^“^^ for a particle with TeV-scale mass. In order to 
make r longer than the age of the universe^, dimension-5 terms (i.e. p = 1) must be effectively 

^ We note that a r comfortably longer than the age of the Universe already also prevents a decaying dark 
matter particle from having an impact on a number of cosmological and astrophysical observations (galaxy and 
cluster formation, type la supernovae, X-ray emissions from clusters, mass-to-light ratios in clusters, cosmic 
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suppressed by A S> Mpi, while dimension -6 operators (i.e. p = 2) are safe for A > 10 ^'^ GeV. 
Therefore in table 1 we also list (in parenthesis) the potentially dangerous dimension-5 opera¬ 
tors. 

One sees that for low n (upper rows of table 1) the multiplets can interact with and decay 
into SM particles in a number of ways. Actually, particles with these quantum numbers already 
appear in a variety of different contexts: e.g. scalar triplets in little-Higgs models; fermion 
or scalar triplet in see-saw models; KK excitations of lepton doublets or of higgses in extra 
dimensional models; higgsinos, sneutrinos, winos in supersymmetric models.^ In all these cases, 
a stable DM candidate can be obtained, but only after suppressing the unwanted decay modes, 
e.g. by invoking some extra symmetry. All the best known WIMP DM candidates need this 
feature (e.g. matter parity in SUSY models [1], T-parity in little-Higgs models [9], KK-parity 
in ‘universal’ extra dimension models [ 10 ]^, etc.). In some cases the overall result of this model 
building activity looks plausible; in other cases it can look more like engineering rather than 
physics. 

On the other hand, known massive stable particles (like the proton) do not decay for a 
simpler reason: decay modes consistent with renormalizability do not exist. DM can be stable 
for the same reason: for sufficiently high n (lower rows of table 1 ), namely 

n > 5 for fermions, and n > 7 for scalars, (2) 

Y is automatically stable because no SM particles have the quantum numbers that allow sizable 
couplings to Y. In other words, DM stability is explained by an ‘accidental symmetry’, like 
proton stability. 

An upper bound on n, 

n < 5 for Majorana fermions and n < 8 for real scalars, (3) 

is obtained by demanding perturbativity of a^^(£'') = ^(Af) — {b2/27i) \nE'/M at E' ~ Mpi, 

where 62 = —19/6 -|- c — l)/36 with c = 1 for fermions and c = 1/4 for scalars. Choosing 
a scale E' much smaller than the Planck mass would not signihcantly relax this upper bound 
on n, in view of the strong rise of 62 with n^, and of the mild logarithmic dependence on E'. 

Therefore table 1 provides a realistically complete list of MDM candidates. 


3 The mass splitting 

At tree level, all the components of X have the same mass M, but loop corrections tend to 
make the charged components slightly heavier than the neutral one. 

microwave background [4, 5]). E.g. measured CMB anisotropies constrain r > 52 Gyr at 95% C.L.. 

^These DM candidates typically behave very differently from the MDM candidates discussed in this paper. 
However, these models have many free parameters, allowing a wide range of different possibilities: in some corner 
of the parameter space they can reproduce the precise phenomenology of the corresponding MDM candidates. 
See e.g. [6, 7, 8] for supersymmetric examples. 

^However, known ways of obtaining chiral fermions from extra dimensions introduce extra structures that 
generically break universality. 
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Actually, in the case of scalar DM there can be an extra source of mass splitting. The 
following extra renormalizable interactions can be present 


^non minimal = ^ - c\h{X*T^X) 


2 


{X*T^Xf 


A' 




Xj 


(4) 


where are SU(2 )l generators in the representation to which R belongs and c = 1 (1/2) for 
a complex (real) scalar. These extra interactions do not induce DM decay (because two X are 
involved, and we assume {X) = 0), and Xx, X'p^ induce no signihcant effects. The coupling A^, 
however, splits the masses of the components of X by an amount 


AM = 


Xhv^\^TI, 


AM 


= Xh ■ 7.6 GeV 


TeV 


(5) 


having inserted {H) = (n, 0) with v = 174 GeV and AT^ = 1. In the following we assume that 
this possible correction to scalar masses is small enough that the mass splitting is determined 
only by loop corrections. We will see that this this is obtained for Xh <0(0.01), which is not 
a unreasonably strong restriction. 

The mass difference induced by loops of SM gauge bosons between two components of X 
with electric charges Q and Q' is explicitly computed as 


Mq — Mo' = 


Q' 


a2M 

Att 


Mz 


(Q" - + iQ- Q')iQ + Q' - 2V) 


M 


/( 


M, 


W-' 


M 


/( 


M 


( 6 ) 


where 


f{r) = 


+r 


2r^In r 
2r^ In r 


2r + (r^ — 4)^/^(r^ + 2) In A /2 for a fermion 
kr + (r^ — 4)^/^ In A /4 for a scalar 


(7) 


with A = (r^ — 2 — ry/r'^ — 4)/2 and sw the sine of the weak angle. The constant k is UV 
divergent, and can be reabsorbed in a renormalization of M and oi Xh ^ Xh + AkYa 2 tan^ 
notice that if V = 0 the UV-divergent term k produces no effect. We will be interested in 
the limit M 3> Mw,z'- for both scalars and fermions the result is well dehned and equal to 

'p _ 

/(r) ~ 2'7ir, such that the mass splitting with respect to the neutral component, Q' = 0, 
becomes: 

Mq - Mo ~ Q{Q + 2Y/ cos 0w) AM for M > Mw, Mz (8) 


where AM is the mass splitting between Q = 1 and Q' = 0 components in the case of zero 
hypercharge, V = 0: 

AM = a 2 Mw sin^ ^ = (166 ± 1) MeV. (9) 

p _^ Q 

As we already observed, when V = 0 the next term in the expansion of /(r) ~ 27rr + {krp 
gives no effect, such that corrections due to hnite M are only of order M^^/M^, and thereby 
practically negligible for the values of M computed in the next section. The values of AM for 
all the relevant MDM candidates are listed in table 1. 

One can easily understand why such loop-induced AM turns out to be not suppressed by 
M (unlike the tree level term of eq. (5)), why the neutral component gets lighter and why spin 
becomes irrelevant for M ^ Mz- The quantum correction AM is dominated by a well-known 
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classical effect: the Coulomb energy [11] stored in the electroweak electric helds. Indeed a 
point-like charge g at rest generates the electric potential ^p{r) = /Airr of one vector 

boson with mass My, such that its Coulomb energy is (in natural units) 


6M 


d^r 




a 

2 


My -l- (UV-divergent constant). 


( 10 ) 


The Coulomb energy on scales M~^ <r<Mz^ gives a common divergent renormalization of 
M. The computable mass difference is the Coulomb energy on scales larger than the SU(2 )l- 
breaking scale r/h> such that the microscopic details of X on much smaller scales ~ M~^ 
become irrelevant. Inserting the W, Z couplings eq. (10) gives AM = a 2 {Mw — c^Mz)/2 which 
is equivalent to eq. (9). 


4 The thermal relic MDM abundance 

Assuming that DM arises as a thermal relic in the Early Universe, we can compute its abundance 
as a function of its mass M. Requiring that X makes all the observed DM, we can therefore 
univocally determine M. 

We can neglect the mass splitting AM M; indeed eq. (25) will conhrm that all DM 
components can be approximatively considered as stable on time-scales comparable to the 
expansion rate at temperature T ~ Mw- The heavier charged components eventually decay 
into the lighter neutral DM, before nucleo-synthesis and giving a negligible entropy release. 
Moreover, following the common practice, we ignore thermal corrections: we estimate that 
they mainly induce thermal mass splittings of order AM ~ {gM)^/M, which likely can also be 
neglected. A more careful investigation of this issue could be worthwhile. 

Thanks to the above approximations we can write a single Boltzmann equation that de¬ 
scribes the evolution of the total abundance of all components A) of the multiplet as a whole. In 
particular, it includes all co-annihilations in the form of J^ij ^A{XiXj —>■ SM particles). Since 
we will hnd that the observed DM abundance is obtained for M^ 3> M|, we can compute 
the relevant thermally-averaged cross-sections in the SU(2)L-symmetric limit. Furthermore as 
usual the freeze-out temperature is Tj ~ M/25 -C M, such that we can keep only the dom¬ 
inant s-wave (co)annihilation processes. The hnal DM abundance can be well approximated 
as [12, 13] 


?7,dm(T) / 180 1 M gxMMpi{aAv) , . 

^(tT “ V WiSMTU- ’ 

where is the number of SM degrees-of-freedom in thermal equilibrium at the freeze-out 
temperature Tf, and s is their total entropy. 

If A is a scalar, annihilations into SM scalars and fermions are p-wave suppressed (notice 
however that the large number of fermions present in the SM could partially compensate this 
suppression), and the dominant annihilation channel is into SU(2 )l U(1)y vector bosons: 
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A'A'* ^ AA. We find^ 


(3 - 4n2 + n^) + 16 + 8glg^yY^(n^ - 1) 

647r M2 gy 


if df is a scalar (12) 


where g^i^ = 2n for a complex scalar and gx = n for a real scalar. 

If df is a fermion, annihilations into gauge bosons are again described by eq. (12), where 
now gx = 4?7, for a Dirac fermion and gx = 2n for a Majorana fermion. Furthermore, now 
also fermion and Higgs hnal states contribute to the s-wave cross section of annihilations plus 
co-annihilations (notice that when H = 0 only co-annihilations are present). The result is 


{(Tav) ~ 


gl {2n^ + 17n^ - 19) + AY^gU^l + 8Y^) + 16g‘lg^Y\n‘^ - 1) 

1287r M2 g;:^- 


if df is a fermion. 


(13) 

Table 1 shows the values of M needed to reproduce all the observed DM.^ Within the standard 
ACDM cosmological model, present data demand = 0.110 ±0.006 i.e. udm/s = (0.40 ± 

0.02) eV/M [15]. In all cases it turns out that M^ M^^, justifying our approximation of 
neglecting SU(2)L-breaking corrections to {(Tav). Since Ddm oc M^, the values of M reported 
in table 1 suffer a 3% experimental uncertainty. We report numerical values indicating only 
the experimental uncertainty because it is difficult to quantify the theoretical uncertainty: it 
can give a comparable contribution, and can be reduced by performing rehned computations. 


In the scalar case, the non-minimal couplings in eq. (4) are generically allowed and give 
extra annihilations into higgses: 


{o' av) extra. — 


|A^|2 + (n2-l)|AH|Vl6 
IOtt M2 gx 


if df is a non-minimal scalar. 


(14) 


(no interference terms are present). The contribution from Xh is negligible for the values that 
allow to neglect its effect on mass splittings (discussed in eq. (5)), and the contribution from 
X'fj is negligible for |A)|^| gy, S'!- We assume that this is the case so that these extra terms are 
negligible for all candidates with electroweak interactions.® In this paper we assume that this is 

"‘We here only make use of 


TrT“T'’ = (5“'’ —±2 - 1), 
12^ 


TrT“T“r^T'’ = — (n2 - l)^ 
16^ ^ 


Trr“T'’r“r'' = — - 5){iY - 1) 


where r“ are the SU(2 )l generators in the representation with dimension n. The neutral DM component in 
the X multiplet has vector-like NC gauge interactions with the Z boson and vector-like CC interactions with 
the boson and with the charged DM^ component. The gauge couplings are 


± g - (It 2±)2 

ffcc ^ • 2 


5NC = 


9Y_ 

Cw 


® As discussed in later works (J. Hisano et ah, hep-ph/0610249 and M. Cirelli et ah, 0706.4071) non pertur¬ 
bative corrections affect the cosmological density in an important way with respect to the perturbative result 
obtained here. 

^However, when A is a neutral scalar singlet, these non-minimal annihilations are the only existing ones: 
the observed amount of DM is obtained for M « 2.2TeV|Ahl (we are assuming M ^ Mz', for generic values of 
M the correlation between M and X'^ was studied in [14]). We stick to the minimal setup so that the singlet 
is not a useful candidate for our purposes. This holds a fortiori for the case of fermionic DM. 
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the case, and here study this issue quantitatively: to compute how much the inferred value of 
M would be affected one simply adds the contribution in eq. (14) to those in eq. (12): the mass 
of the low n candidates are more affected than the high n ones, because the gauge-mediated 
annihilation cross sections are relatively less important for the formers. For example, a large 
X'jj = 1 would increase the predicted value of M by a factor 2.4 for n = 2, by 20% for n = 3, 
by 2% for n = 5, and by 0.5% for n = 7. Note that Xh ~ 92 and X'^ ~ gy are the values 
predicted by some solutions to the hierarchy problem, such as supersymmetry and gauge/higgs 
unihcation. 

A smaller M is allowed if X gives only a fraction of the total DM density. In the more general 
situation where many MDM multiplets are present, to a good approximation their abundances 
evolve independently, such that the observed DM abundance is reproduced for lower values 
of DM masses. In this more general situation, the M values in table 1 must be reinterpreted 
as upper bounds on M. For example, with N identical families of the same X particle their 
common mass becomes times lower. 

Finally, we notice that the combination of a precise collider measurement of M with a precise 
measurement and computation of the DM abundance, could allow to test the expansion rate at 
a freeze-out temperature around the electroweak scale. This is a particularly interesting region: 
one can test if the electroweak vacuum energy behaves just like as a fine-tuned constant, or if 
some sort of relaxation mechanism of the cosmological constant takes place. 

5 Low-energy MDM signals 

Direct detection 

Searches for elastic DM collision on nuclei A/” provide strong constraints. MDM candidates with 
D % 0 have vector-like interactions with the Z boson that produce spin-independent elastic 
cross sections 

jYf2 

a(DMA/' ^ DMAA) = c ^ Y\N - (1 - 4s^)Z)2 (15) 

ZTT 

where c = 1 for fermionic DM and c = 4 for scalar DM [16]; Z and N are the number of 
protons and of neutrons in the target nucleus with mass M^/-, We are assuming M ^ M/y. 
For all DM candidates with D % 0, this elastic cross section is 2 3 orders of magnitude 

above present bounds [17]. Such MDM candidates are therefore excluded, unless minimality 
is abandoned in appropriate ways that allow to avoid the experimental limit [6, 7]; this for 
instance happens is the well-known case of the Higgsino: the mixing with Majorana gauginos 
can split its components by an amount 5m such that the lightest one becomes a Majorana 
fermion which cannot have a vector-like coupling to the Z boson. In our case, if a 5m larger 
than the DM kinetic energy and smaller than Mq—Mq is generated by some similar non-minimal 
mechanism, it would affect only the direct detection DM signals such that all candidates with 
D % 0 become still allowed: therefore they are included in table 1. However, in the following 
we insist on minimality and focus mainly on the candidates with Y = 0. 

MDM candidates with D = 0 have vanishing DMA/” and DM DM cross sections at tree level, 
but are accompanied by a charged component heavier by AM = 166MeV (see eq. (9)), that 
produces various possible effects. 





Figure 1: One loop DM/quark scattering for fermionic MDM with F = 0. Two extra graphs 
involving the four particle vertex exist in the case of scalar MDM. 


An elastic cross section on nuclei is generated at loop level via the diagrams in fig. 1. An 
explicit computation of the relevant one-loop diagrams is needed to understand qualitatively and 
quantitatively the result. We find that non-relativistic MDM/quark interactions of fermionic 
MDM with mass M S> M\y S> are described by the effective Lagrangian 




(n 


(l±2y)^ 


nan 


16My[r 


E 


<? L 


'M^ 


+ 


m] 


r)[AA]mg[gg] 


2 




(16) 


where the -|- (—) sign holds for down-type (up-type) quarks q = {u, d, s, c, b, t}, nih is the Higgs 
mass and nig are the quark masses. The first operator gives dominant spin-independent effects 
and is not suppressed by M; the second operator is suppressed by one power of M and gives 
spin-dependent effects. Parameterizing the nucleonic matrix element as 


(A^l ^ mqqq\N) = fniN (17) 

<? 

where mpf is the nucleon mass, the spin-independent DM cross section on a target nucleus A/" 
with mass Mj^ is given by 

...(DMA^ ^ DMA^) = (18) 

In the case of scalar MDM we find in the relevant non-relativistic limit: an M-independent 
contribution to asi equal to the fermionic result of eq. (18); an UV-divergent effect suppressed 
by M that corresponds to a renormalization of lApliPp operators (that can produce a much 
larger usi if present at tree level); no spin-dependent effect. 

Assuming nih = 115 GeV and / ~ 1/3^ (QCD uncertainties induce one order of magnitude 
uncertainty on usi) we plot in fig. 2 the MDM prediction for the standard nucleonic [1, 19, 22] 

^To properly compute nuclear matrix elements one must keep quarks off-shell, finding several operators that 
become equivalent on-shell [19]: 

mq[XX][qqf [XX][qi^q], ■^[Xid^'y^X][qi{dg"f^ + -^^)q],.. . 

Summing over all quarks the matrix elements are / « (0.3 -:- 0.6) [20] for the first operator, / k, 1.2 for the 
third operator, while the matrix element of the second operator is unknown. In our computation, only the 
first operator contributes to the SI effects suppressed by the higgs mass, while the other SI effects arise from 
a combination of the various operators in proportion 0 : — I : 2. Therefore cancellations are possible. We do 
not fully agree with result of a previous computation [21], performed for the fermionic supersymmetric DM 
candidates: wino (n = 3, Y = 0) and Higgsino (n = 2, |F| = 1/2). 
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DM mass in GeV 


Figure 2: Spin-independent cross sections per nucleon of MDM candidates assuming 
the matrix element f = 1/3. Fermionic MDM candidates are denoted as blue circles; scalars 
as red diamonds. The fully successful candidates (automatically stable Y = 0 multiplets with 
n = 5,7) are represented by filled symbols. We also show, as large empty symbols, the Y = 0 
candidates that reguire a stabilization mechanism and, with smaller symbols, those with Y ^ 0 
(viable if an additional mechanism forbids a much larger cross section mediated by a Z boson). 
The dashed lines indicate the sensitivity of some future experiments [18]. The cloud indicates 
the range of values favoured by a minimal SUSY model. 


cross section parameter, showing that it requires experimental sensitivities about 3 orders of 
magnitude larger than the current ones, but within the sensitivity of experiments currently 
under study. The annual modulation effect of the DAMA/Nal experiment [23] cannot be 
explained by MDM candidates, since they have too large masses and too small cross sections 
with respect to the properties of a WIMP compatible with the effect. We also show the signal 
rates for the MDM candidates with Y ^ 0, assuming that they are resurrected by splitting their 
neutral components. In this case, the W contribution of eq. (16) to the effective Lagrangian 
relevant for direct DM signals becomes 2 times lower, and there is an extra Z contribution 
given on-shell by 


c^Z _ 


C^Mz q 


3 [A A] 


mq[qq] 


{dLq + 9Rq) 

2M| 


+ 


4m? 


M 


(glq + 9lq) 


where 9l = Ts — and gn = —Qs"^ describe Z/quark couplings. 

For comparison, the cloud in £g. 2 represents a typical range of values found in scannings over 
the MSSM parameter space [24], restricted by present bounds and by naturalness assumptions. 
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CC direct detection? 


MDM could be probed via the tree level CC process DMA? —> A/'^DM^. In ordinary conditions 
this process is kinematically forbidden, because neither the DM kinetic energy, nor the A? 
kinetic energy, nor the nuclear mass difference is large enough provide an energy larger than 
the AM of eq. (9). In these conditions a virtual DM^ can be produced slightly off-shell, but 
the cross section for the resulting multi-body process DMA? —> DMA?^e^^I7^ is negligibly 
small. Production of on-shell DM^ could take place in environments, like cluster of galaxies or 
supernovae, where DM or nuclei have enough kinetic energy; but we do not see how this can 
provide a detectable signal. 

We here explore the possibility of accelerating nuclei in a controlled experiment: a dedicated 
discussion seems worthwhile because prospects of searching for this unusual signal seem not 
much beyond what can be considered as realistic. Protons p or nuclei A? can be accelerated 
up to energy Ep > AM (while collider experiments need a much larger energy Ep > M) and 
stored in an accumulator ring. The DM background acts as a diffuse target such that there 
is no need for focusing the beam. The cross section is more easily computed in the limit 
AM, rUp Ep M. Averaging over initial colors and polarizations we hnd the following 
partonic cross section that applies for both fermionic and scalar MDM: 


(j(aDM —>■ a^DM^) = cxo 


— 1 


ln(l + 4EVM^) 


(To = 


TT 


= 1.110“^^cm^ 


(19) 


where a is any quark or anti-quark (or electron) with energy E. Therefore the pDM —> n DM"*" 
cross section approaches the constant value 3 (To(u^ — l)/4 for Mw ^Ep<M. Increasing the 
proton energy Ep above Mw increases the DM^ velocity, f3 ~ max(10“^, {Ep/My^'^), but does 
not increase the cross section. Such a high cross section, not suppressed by powers of M, arises 
because this is a non-abelian Coulomb-like process. For comparison, SUSY DM candidates have 
comparable cross sections (lower if the neutralino is dominantly bino) and chargino/neutralino 
mass splittings AM possibly in the 10 GeV range.® 

At energies Ep -C Mw the cross section decreases as {Ep/MwY, and when Ep<mp one 
should switch from a partonic to a nucleonic description. Nuclear cross sections are not coher¬ 
ently enhanced for the AM values suggested by MDM. The event rate is 


dN _ ^ pdm _ 10 Np pdm TeV a 

dt ^ M year 10^° 0.3 GeV/cm® M 3 (To 


( 20 ) 


where pdm is the local DM density and e is the detection efficiency, related e.g. to the fraction 
of beam that can be monitored. A circular accumulator ring that employs a magnetic field B 
has radius R = Ep/qpB = 21 m(i?p/10 GeV)(10 Tesla/i?). Proton drivers currently planned for 
neutrino beam experiments can produce more than 10^® protons per second, and accumulating 
Np ~ 10^° protons is considered as possible. 

®A better experiment might be possible in the case of supersymmetric neutralino DM, but only after that 
sparticle masses will be precisely measured: an electron beam can be stored at the energy that makes the 
process electron neutralino ^ selectron resonant, with the production of electrons of known energy [25]. In this 
respect, we also note that plans for neutrino factories are independently proposing muon beams (less plagued 
by synchrotron losses than electron beams) that could have the desired high intensity (10^^ muons stored per 
year?), time-scale (2020?) and energy (few tens of GeV?) for a muon neutralino —> smuon experiment; for 
detection, a muon signal is better than an electron signal. 
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The main problem seems disentangling the signal from the beam-related backgrounds. This 
issue crucially depends on how DM^ behaves; e.g. an enough long lived DM^ decaying into 
would give a clean enough signal. In section 6 we compute the DM"*" life-time, hnding that it 
is long lived, but not enough. 

Indirect detection 

DM DM annihilations can have larger cross-sections. Scalar DM annihilations into higgses can 
arise at tree level. Furthermore, as noticed in the context of SUSY models with multi-TeV LSP 
mass® annihilations into gauge bosons of DM particles with non relativistic velocity (3 < Oem (in 
our galaxy (3 ~ 10“^) can be strongly enhanced by the presence of quasi-degenerate charged 
SU(2 )l partners [26]. For MDM candidates with Y = 0, by performing a naive perturbative 
computation of the relevant Feynman diagrams (see [27] for analogous computations in the 
context of supersymmetry), we hnd 

a{DMDM^W+W-)v = (n^ - ( 21 a) 

at tree level, and 

TTcy^ 

a(DMDM^ 77 )u = (n^ - l)^-^ (21b) 

at one loop. Analogous results hold for Z 7 and ZZ final states. However, the DM®DM® 
system can happen to have mass equal to the DM’'"DM“ or DM’'"^DM systems: the mass 
difference AM ~ a 2 Mw can be compensated by the binding energy Eb ~ a|M of the two- 
body state. This happens for specific values M* ~ My//a 2 of M: if M ^ M* the cross sections 
cr(DMDM —> AA) get enhanced by a factor 0{1 — M/M*)“^ and acquire comparable values 
for all vector bosons A [26]. Unlike in the supersymmetric case, in the MDM case AM and 
consequently M* are univocally predicted. We find resonances at^® 


n 


M, 

in TeV 

3 

2.5 

9.8 


5 

1.8 

3.3 

6.6 ... 

7 

.74 

1.6 

2.9 3.7 ... 


( 22 ) 


®Iii this corner of its parameter space SUSY does not naturally provide electroweak-symmetry breaking: 
cancellations among individual corrections to the squared Higgs mass which are thousand of times larger than 
their sum are needed. In this region DM DM annihilations are negligibly affected by all other sparticles outside 
the DM multiplet, restricted to be higgsino (n = 2) or wino (n = 3). 

^°We here outline the computation, performed along the lines of [28, 26]. Using the case n = 5 as an 
example, the Hamiltonian that describes the non-relativistic canonically normalized {DM''''*'DM , DM'*'DM“, 

DM®DM®} two-body states is 


H = - 


Yl 

Jd 


AMdiag(8,2,0) - 


, I dttem + 4a2Cwe 
2a2e-^'^^ 
0 


2026 "“'^’' 
Uem + 

3U2a2e-^^«'’' 


0 

3V2a2e-^'^^ 
0 


One recognizes the kinetic energy in the center-of-mass frame, the contributions due to the mass difference AM, 
the Coulomb energy and its non-abelian generalizations. We numerically found the values of M that 

give rise to s-wave eigenstates with H = 0, relevant for DM°DM® annihilations. 
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In the case n = 3 the hrst resonances happens to be close to the value of M suggested by the 
DM abundance, giving a (9(100) enhancement. 

DM DM annihilations can produce different types of signals relevant for indirect DM 
searches. We do not perform more precise computations of the relevant particle-physics because 
the predicted rates are signihcantly affected by astrophysical uncertainties concerning the DM 
distribution in our Galaxy. 

When DM annihilations occurs in the inner core of the Earth or the Sun, producing a 
neutrino flux to be searched at neutrino telescopes, the signal is strictly linked to the situation 
which occurs for direct detection. In this case the signal depends on the scattering cross section 
on the nuclei of the capturing body (Earth or Sun). We estimate that the MDM candidates 
compatible with direct detection (plotted in hg. 2) provide a neutrino signal which is about 
two orders of magnitude lower than current sensitivities [29]. Large area detectors, like e.g. 
ICECUBE [30], might access the upgoing muon signal for these DM candidates. The neutrino 
energy and flavour spectra can be computed in terms of particle physics along the lines of [31]. 

In the case of DM annihilation inside the Galactic halo, the signals consist in fluxes of 
antiprotons, antideuterons, positrons and gamma-rays. We estimate that an annihilation cross 
section around 10“^^ cnU/s is needed in order to provide an antiproton signal which can emerge 
over the expected background for antiproton kinetic energies lager than about 100 GeV, and 
without falling in conflict with the lower energy data [32]. The off-resonance cross section of 
eq. (21a) can only reach cm^/s for n = {5,7}. Signals at the level of the expected 

background can be approached for DM masses that lie within ~ 10% of the resonant values 
M* discussed above. In this case, forthcoming experiments like PAMELA [33] and AMS [34] 
have a chance to observe a signal. However, we recall that the antiproton signal is affected 
by a large astrophysical uncertainty, of the order of a factor of a few [32]. This may induce a 
more favourable situation for the antiproton signal, especially in the case of a thick conhnement 
region in the Galaxy. As for positrons, again enhanced cross sections for M close to the resonant 
values M* can provide signals at the level of detectability [26]. 

A cross section at the level of providing a signal in the antiproton channel would also give a 
antideuteron signal [36] comparable to its background for energies larger than about 50 GeV. 
However, planned experiments like GAPS [35] are sensitive to the more relevant low-energy 
tail [36]; DM annihilations into gauge bosons are not a favourable channel for antimatter 
production [32, 37]. Therefore, the possibility to detect antideuterons from MDM looks rather 
difficult. 

Finally, in the case of gamma-rays, we extimate that a cross section at the level of 
cm^/s can allow detection of a signal from the galactic center for the GLAST experiment [38, 39], 
in the energy range around 100 GeV and assuming a Navarro-Frenk-White density prohle. Gross 
sections one order of magnitude smaller can be accessed by future water Cerenkov detectors 
like VERITAS [40, 39]. Again, also in this case only a resonantly enhanced cross section can 
provide a detectable signal; in such a case a signal coming from external Galaxies like M31 
could be accessed from experiments like EGRET and VERITAS [39]. 
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6 High energy MDM signals 


MDM multiplets give one loop ‘universal’ corrections to electroweak precision data affecting 
only the W, V parameters [41] as 


W 


02 -1 
^ ^^6071 12 


Y 


c 


dOvT 


(23) 


where c = 1 for fermionic DM and c = 1/4 for scalar DM. For the masses M computed in 
section 4 one gets a negligibly small IF, F ~ 10“^. 

Direct pair production of MDM particles is a more promising future signal. We focus on the 
LHC collider, that will collide pp at ^/s = 14 TeV. For F = 0 the partonic total cross sections 
(averaged over initial colors and spins) for producing all DM components are 


^ud ^du 2(J 


uu 


2a 


dd 


-1) 

13824 TTS 


j3‘^ if F is a scalar 
3-/3^ if F is a fermion 


(24) 


where the subscripts denote the colliding partons, and /? = y 1 — AM'^/s is the DM velocity 
with respect to the partonic center of mass frame. Production of non-relativistic scalars is p- 
wave suppressed in the usual way. In table 1 we show the number of pp —*■ DMj DMj X events 
(here X denotes any other particles). We computed cross sections including the hypercharge 
contribution, summed over all DM components, and assumed an integrated luminosity equal to 
100/fb, that each detector will accumulate in 10^ sec (i.e. one collider-year) if LHC will deliver 
its planned luminosity, C ~ 10^°/cm^sec. The event rate depends signihcantly on the precise 
value of the DM mass M, but surely some MDM candidates are too heavy for LHC. A collider 
with beam energy 2-^4 times higher could fully test all MDM candidates. Possible upgrades 
of LHC luminosity and magnets are discussed in [42]. 

MDM has a clean signature, that allows discovery even if only a few events are seen. The 
small mass splitting among the DM components makes too hard to tag the missing energy 
carried away by neutral DM particles, but also makes charged MDM component (s) enough 
long-lived that they manifest in the detector as charged tracks. Irrespectively of the DM spin 
the life-time of DM^ particles with F = 0 and n = {3, 5, 7,...} is r ~ 44 cm/(n^ — 1) and the 
decay channels are 


DM^ ^ DM°7r± 
DM^ ^ 

DM^ ^ DMV^H 


r. = (n^ - 1) 


GlV^AAPfi 


r, = (n- 


4:71 


' 1 - 


mt 


AM2’ 


bOvr^ 


Tf, = 0.12 P, 


BR^ = 97.7% 

BRe = 2.05% 
BR^ = 0.25% 


(25) 

having used the normalization = 131 MeV [43] and the AM of eq. (9), which accidentally 
happens to be the value that maximizes BR,r- The DM"*" life-time is long enough that decays 
can happen inside the detector. On the contrary, the faster decays of DM^^ particles (present 
for n > 5) mostly happen within the non-instrumented region with few cm size around the 
collision region. Measurements of r and of the energy of secondary soft pions, electrons and 
muons constitute tests of the model, as these observables negligibly depend on the DM mass M. 
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On the contrary, measurements of the total number of events or of the DM velocity distribution 
would allow to infer its mass M and its spin. Although SM backgrounds do not fully mimic 
the well defined MDM signal, at an hadron collider (such as the LHC) their rate is so high that 
some experimentalists consider impossible triggering on the MDM signal. 

Notice that extra SU(2 )l multiplets that couple (almost) only through gauge interactions 
tend to give a LHC phenomenology similar to the one discussed above, irrespectively of their 
possible relevance for the MDM problem. On the contrary, DM candidates like neutralinos are 
often dominantly produced through gluino decays, such that DM is accompanied by energetic 
jets rather than by charged tracks. 

7 Conclusions 

We extended the Standard Model by adding a spin-0 or spin-1/2 n-tuplet of SU(2)l with 
hypercharge Y that only has gauge interactions and mass M. Some multiplets contain neutral 
components, that are potential Dark Matter (DM) candidates. 

• Multiplets with Y ^ 0 are already excluded by direct DM searches. They can be res¬ 
urrected by introducing non-minimal mechanisms that prevent Z-mediated DM/nuclei 
coupling, e.g. by appropriately mixing their neutral components with a singlet. 

• Multiplets with Y = 0 and odd n = {3, 5, 7,...} contain allowed DM candidates. 

— For n = 3 one needs to impose DM stability by hand. 

— For n > 5 the stability is instead automatically guaranteed by renormalizability, 
much alike proton stability. 

The set of interesting candidates is bounded by n< 7 in order to avoid Landau poles in 02 - 

Gauge interactions are spontaneously broken and thereby induce a non-trivial and peculiar 
Minimal DM (MDM) phenomenology, fully computable in terms of a single unknown parameter: 
the DM mass M. Electroweak breaking effects induce a mass splitting AM ~ aMw among 
the components of any given multiplet, making the neutral component lighter than the charged 
components. Assuming that only one MDM multiplet is present, its mass M is determined 
by the request that its relic thermal abundance equals the observed DM abundance. Co¬ 
annihilations play a crucial role, giving M ~ few TeV. Since M 3> My/-, various MDM properties 
depend dominantly only on the MDM gauge charge, while the microscopic MDM properties 
(such as their spin) become irrelevant. 

The simplest fully successful MDM candidate is a fermionic SU(2 )l quintuplet with mass 
M 4.4 TeV. MDM candidates are listed in table 1: some are fully successful (automatically 
stable and consistent with DM searches), others require a stabilization mechanism (e.g. the 
wino-like candidate) or a way to elude the bounds from direct DM searches (we list only those 
with n < 4), or both (e.g. the Higgsino-like candidate). If multiple MDM multiplets exist, all 
their masses become lighter than in table 1: e.g. 42% lighter in presence of 3 identical families 
of a single multiplet, signihcantly increasing the number of events expected at LHC. 

MDM multiplets contain charged components, slightly heavier than the neutral DM com¬ 
ponent. For V = 0 the charged DM^ is AM = 166 MeV heavier and has a life-time r = 


15 


44cm/(n^ ~ 1); giving a clean displaced-vertex signature at colliders. The branching ra¬ 
tios are predicted to be BR(DM± ^ 7r±DM°) = 97.7%, BR(DM± ^ DM°e±¥e) = 2.05%, 
BR(DM^ —= 0.25%. We computed the event rate at LHC, finding that LHC 
cannot probe all MDM multiplets, being more sensitive to the ones with lower n (and, if mul¬ 
tiple multiplets are present, to the ones that would give subdominant contributions to the DM 
density). 

On the contrary direct DM searches are more sensitive to higher n MDM multiplets (and to 
ones that dominate the DM density). Indeed one-loop diagrams generate a spin-independent 
MDM/nucleus cross section parameter usi ~ 10“^^(?7,/5)'^ cm^ (up to a QCD uncertainty of 
about one order of magnitude). As illustrated in fig. 2, this is within the sensitivity of future 
experiments. 

DM DM annihilations in the galactic halo into vector bosons can be resonantly enhanced, 
giving indirect DM signals, only for DM masses close to the values listed in eq. (22). CC 
production of DM^ occurs at tree level with a much larger cross section, that can exceed 
10“^'^ cm^, but in ordinary situations it is forbidden kinematically. We also discussed prospects 
for attempting CC direct DM detection by accelerating an intense nuclear beam. 
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